A proposal for testing Quantum Gravity in the lab 



O 
(N 

< 



Oh' 
CD 



> 

rn 

O 



X 



Ahmed Farag Ali ^,* Saurya Das ^,^ and Elias C. Vagenas '^^ 

^ Theoretical Physics Group, Department of Physics, 

University of Lethbridge, 4401 University Drive, 

Lethhndge, Alberta, Canada TIK 3M4 

and 

^ Research Center for Astronomy and Applied Mathematics, 

Academy of Athens, 

Soranou Efessiou 4, GR-11527, Athens, Creece 

Attempts to formulate a quantum theory of gravitation are collectively known as quantum gravity. 
Various approaches to quantum gravity such as string theory and loop quantum gravity, as well as 
black hole physics and doubly special relativity theories predict a minimum measurable length, or a 
maximum observable momentum, and related modifications of the Heisenberg Uncertainty Principle 
to a so-called generalized uncertainty principle (GUP). We have proposed a GUP consistent with 
string theory, black hole physics and doubly special relativity theories and have showed that this 
modifies all quantum mechanical Hamiltonians. When applied to an elementary particle, it suggests 
that the space that confines it must be quantized, and in fact that all measurable lengths are quan- 
tized in units of a fundamental length (which can be the Planck length) . On the one hand, this may 
signal the breakdown of the spacetime continuum picture near that scale, and on the other hand, 
it can predict an upper bound on the quantum gravity parameter in the GUP, from current obser- 
vations. Furthermore, such fundamental discreteness of space may have observable consequences at 
length scales much larger than the Planck scale. Because this infiuences all the quantum Hamilto- 
nians in an universal way, it predicts quantum gravity corrections to various quantum phenomena. 
Therefore, in the present work we compute these corrections to the Lamb shift, simple harmonic 
oscillator. Landau levels, and the tunneling current in a scanning tunneling microscope. 



I. INTRODUCTION 



An intriguing prediction of various theories of quantum 
gravity (such as string theory) and black hole physics is 
the existence of a minimum measurable length. This has 
given rise to the so-called generalized uncertainty prin- 
ciple (GUP) or equivalently, modified commutation re- 
lations between position coordinates and momenta. The 
recently proposed doubly special relativity (DSR) theories 
on the other hand, also suggest a similar modification of 
commutators. The commutators that are consistent with 
string theory, black holes physics, DSR, and which ensure 
[xi,Xj\ = = [PitPj] (via the Jacobi identity) have the 
following form [1] (see Appendix) ^ 



[xi,pj\ = ihi 5ij-a{p5i.j + ^-^ j + a'^{p'^5^j + 3piPj)j(l) 



where a = ao/Mpic — ao£pi/h, Mpi — Planck mass, 
ipi « 10~^^ m = Planck length, and Mpic^ = Planck 
energy « 10^^ GeV . 



In one dimension, Eq.(l) gives to Oio?) 

AxAp > -[l-2a <p> +4a^ < p^ >] 



> 



+ 4:a^]Ap^+4a^{p)^-2ay^Jj^ 



(2) 



Commutators and inequalities similar to (1) and (2) were 
proposed and derived respectively in [2-9]. These in 
turn imply a minimum measurable length and a maxi- 
mum measurable momentum - the latter following from 
the assumption that Ap characterizes the maximum mo- 
mentum of a particle as well [10], and also from the fact 
that DSR predicts such an maximum (to the best of our 
knowledge, (1) and (2) are the only forms which imply 
both) 



Aa; > (Aa:) 
Ap < (Ap) 



mm ~ (XqIpi 

Mpic 



ao 



Next, defining (see Appendix) 



Xi Xq2 



Poi (1 - apo + 2a Pq) 



(3) 

(4) 

(5) 
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^ The results of this article do not depend on this particular form of 
GUP chosen, and continue to hold for a large class of variants, so 
long as an 0{a) term is present in the right-hand side of Eq.(l). 



with XqijPqj satisfying the canonical commutation rela- 
tions [xoijPoj] = ifi 6ij, it can be shown that Eq.(l) is 
satisfied. Here, poi can be interpreted as the momentum 
at low energies (having the standard representation in 
position space, i.e. poi — —ihd/dxQi) and pi as that at 
higher energies. 



It is normally assumed that the dimensionless parameter 
Uq is of the order of unity, in which case the a depen- 
dent terms are important only when energies (momenta) 
are comparable to the Planck energy (momentum), and 
lengths are comparable to the Planck length. However, if 
we do not impose this condition a priori, then this may 
signal the existence of a new physical length scale of the 
order of ah = aotpi. Evidently, such an intermediate 
length scale cannot exceed the electroweak length scale 
~ 10"'^^ £p; (as otherwise it would have been observed) 
and this implies that ag ^ 10^^. 
Using Eq.(5), a Hamiltonian of the form 



H = |- + n^i 



can be written as 



H 

where H^ 

and Hi 



Ha + Hi + Oia^) 

f- + V{r) 
2m 



Po 



2m 



Po 



(6) 

(7) 
(8) 

(9) 



Thus, we see that any system with a well-defined quan- 
tum (or even classical) Hamiltonian Hq is perturbed by 
Hi, defined above, near the Planck scale. Such correc- 
tions extend to relativistic systems as well [11], and given 
the robust nature of GUP, will continue to play a role ir- 
respective of what other quantum gravity corrections one 
may consider. In other words, they are in some sense uni- 
versal. 

The relativistic Dirac equation is modified in a similar 
way and confirms the main results of our paper [11]. In 
this paper, we first study the effects of the above GUP- 
corrccted Hamiltonian to a particle in a box, to 0{a) in 
Sec. lA, and to 0(0^) in Sec. IB, and show that they 
lead to virtually identical conclusions. In Sec. II, we 
study the effects of GUP-corrected Hamiltonian to the 
Landau levels. In Sec. HI, we calculate the corrections 
due to GUP in the context of a simple harmonic oscilla- 
tor. In Sec. IV, we study the effects of GUP on the Lamb 
shift. Furthermore, we compute the GUP corrections on 
the tunneling current in a scanning tunneling microscope 
for a step potential in Sec. V and for a potential barrier 
in Sec. VI. Finally, we summarize our results in the con- 
cluding section. 



A. Solution to order a 

In this subsection, we briefly review our work in [1] . The 
wave function of the particle satisfles the following GUP- 
corrected Schrodinger equation inside the box of length 
L (with boundaries at x = and x ~ L), where V{r) = 
(outside, U = oo and V' = 0) 



which is now written, to order a, as 

(fi) + k^i) + 2iahStP = 



(11) 



where d" stands for d^/dx" and k ~ ^2mE/h?. A trial 
solution of the form ijj — e™^ yields 

m"^ + P + 2iahm^ = (12) 

with the following solution set to leading order in a: m — 
{ik' , -ik" ,i/2ah}, where k' = k{l + kah) and fc" = fc(l- 
kah). Thus, the general wavefunction to leading order in 
ipi and a is of the form 



iP ^ Ae'*^'^ + Se-**^""^ + Ce*="/2"fi. 



(13) 



Although the flrst two terms can be considered as per- 
turbative corrections over the standard solutions, the 
appearance of the new oscillatory third term is note- 
worthy here, with characteristic wavelength Anah and 
momentum l/4a = Mpic/Aao [which is Planckian for 
ao — 0{1)]. This can be termed a nonperturbative so- 
lution as the exponent contains 1/a and results in the 
new quantization mentioned above. Note that however, 
as explained in [1] and [11], C scales as a power of a, and 
the new solution disappears in the a — > limit. 
Imposing the appropriate boundary conditions, i.e. ■0 = 
at X — 0, L, with A assumed real without loss of gen- 
erality, we get, to leading order, the following two series 
of solutions (C= ]C]e-'^c^): 

TT + 2qn + 29c ^ Ptt + 29c (14) 



2ah 

L 
2ah 
P 



2aoipi 
L 

2ao£pi 
2q±neN. 



2qTT = pi: 



(15) 



Hi} = Ei} 



(10) 



These show that there cannot be even be a single particle 
in the box, unless its length is quantized as above. For 
other lengths, there is no way to probe or measure the 
box, even if it exists. Hence, effectively all measurable 
lengths are quantized in units of ao^p;. We interpret this 
as space essentially having a discrete nature. Note that 
the above conclusion holds for any unknown but fixed 
6c, which, however, determines the minimum measur- 
able length, if any. It is hoped that additional physically 
motivated or consistency conditions will eventually allow 
one to either determine or at least put reasonable bounds 
on it. 

The minimum length is « aoipi in each case. Once 
again, if ao ~ 1, this fundamental unit is the Planck 
length. However, current experiments do not rule out 
discreteness smaller than about a thousandth of a Fermi, 
thus predicting the previously mentioned bound on ag- 
Note that similar quantization of length was shown in 
the context of loop quantum gravity in [12]. 



B. Solution to order a 

We extend the previous solution to include the a^ term 
in one dimension. Working to 0(0^), the magnitude of 



the momentuni at high energies as given by Eq.(5) reads 

p = po{l - apo + 2a^pl) . (16) 

The wavcfunction satisfies the following GUP-corrected 
Schrodinger equation 

dV + k'^i' + 2ihad^tP - hh^a^d^iP = (17) 



where k = ^2mE/h? and d" = rf"/dx". 
Substituting i/'(x) = e™^, we obtain 

m^ + fc2 + 2iahm^ - h{aKfm^ = (18) 

with the following solution set to leading order in o?: 
m = {ik', -ik", |i|, ^2±i}, where k' = k{l + kah) and 
k" = k{l — kah). Thus, the most general solution to 
leading order in £p; and a^ is of the form 

+ i:)e(-2+i)^/5"'i . (19) 

Note again the appearance of new oscillatory terms, with 
characteristic wavelength lOnah, which as before, by 
virtue of C and D scaling as a power of a, disappear 
in the a — > limit. In addition, we absorb any phase of 
A in ip so as A to be real. The boundary condition 

V'(O) = (20) 

implies 

A + B + C + D = (21) 

and hence the general solution given in Eq.(19) becomes 



ipix) = 2iAsm{kx)e'°'''''''' -{C + D)e-'''"'' 



(22) 



If we now combine Eq.(22) and the remaining boundary 
condition 



^{L) =. 
we get 

2iAsm{kL) = ((7 + £))e-»["fe'sL+fe"L] 



(23) 



Ce'soR + De Sen 



eihe-''^'''''^ .(24) 



We can consider the exponentials e^*"*^ ''■^ « 1, oth- 
erwise, since they are multiplied with C or D, terms 
of higher order in a will appear. Therefore, we have 



-ikL 



2iAsm{kL) ^ IC'le^'"^' + \D\e-"'" 



h p 5a?i 



.(25) 



Now, equating the real parts of Eq.(25) (remembering 
that A e M), we have 



= \C\ cos{0c + kL) + \D\ cos{0D + kL) 
— gs^ \C\ cos{0c — -^^) — e^^fi" \D\ cos(0£) 



L 



(26) 



5ah bah 

Note that the third term in the right hand side dominates 
over the other terms in the limit a — > 0. Thus we arrive 
at the following equation to leading order 



cos{L/bah— 6c) = 



(27) 



This implies the quantization of the space by the follow- 
ing equation 



L 

5ah 



(2p+l)-+9c , pe N 



(28) 



Once again, even though the a^ term has been included, 
the space quantization given in Eq.(28) suggests that the 
dimension of the box, and hence all measurable lengths 
are quantized in units of aoipi, and if ag ~ 1, this fun- 
damental unit is of the order of Planck length. And as 
before, the yet undetermined constant 6c determines the 
minimum measurable length. 



II. THE LANDAU LEVELS 

Consider a particle of mass m and charge e in a constant 
magnetic field B ~ Bz, described by the vector potential 
A = Bxy and the Hamiltonian 



2m 2m mi 



2h2 



e'B 
2m 



(29) 
(30) 



Since poy commutes with H , replacing it with its eigen- 
value Sfc, we get 



^n 



^ + Imujl I X 

2m 2 \ moj. 



hk 



(31) 



where ujc ~ eB/m, is the cyclotron frequency. This is 
nothing but the Hamiltonian of a harmonic oscillator 
in the x direction, with its equilibrium position given 
by xq = hk/mtOc- Consequently, the eigenfunctions and 
eigenvalues are given, respectively, by 

V'fe,„(a;,2/) = e'''y4>n{x~xo) (32) 

En = hujJn+-] , neN (33) 



where (/)„ are the harmonic oscillator wavefunctions. 
The GUP-corrected Hamiltonian assumes the form [9] 



H = 



^{p'-'^y-^{p'-''^y 



5a^ 
2m 



Po 



eA 



= Hq — V 8™ a Hq 



10 a^ m iJo^ (34) 



where in the last step we have used Eq.(29). Evi- 
dently, the eigenfunctions remain unchanged. However, 
the eigenvalues are shifted by 



AR 



n{GUP) 



'%ra uHq + 6 a tuHq |(/)„) 



'8m a (hujc) ^ [ 1^+ t: 



+ 10 ma^ihujcf [n + 



(35) 



which can be written as 



A£;, 



n{GUP) 



(0) 



— V 8™ a{hjJc) ^ \n+ - 



+ 10 ma^ihujc) n + 



E, 



For n=l, we obtain the following relation 



^-£'l(GC/P) 



El 



(0) 



+ 



\/l2m (hujc)^ 
Mpic 

15 m {fvjjc) 2 



"0 



(37) 



For an electron in a magnetic field of lOT, cj^ ~ IQ^GHz 



AEi 



(GUP) 



E- 



(0) 



-10-2% 



10-^^al 



(38) 



Thus, quantum gravity/GUP does affect the Landau lev- 
els. However, once again, assuming ao ^ 1 renders the 
correction too small to be measured. Without this as- 
sumption, due to an accuracy of one part in 10^ in direct 
measurements of Landau levels using a scanning tunnel 
microscope (STM) (which is somewhat optimistic) [13], 
the upper bound on ao becomes 



ao < 10 



23 



(39) 



Note that this is more stringent than the one derived in 
previous works [9]. 



Employing time-independent perturbation theory, the 
eigenvalues are shifted up to the first order of a by 



A^, 



GUP 



= {^n\Hi\^n) 



(41) 



where V'n a-re the eigenfunctions of the simple harmonic 
oscillator and are given by 



where 



2"n! 



nh 



-H^J^x) (42) 



2 rf" 

Hn{x) = (-l)"e" -^e' 



(43) 



are the Hcrmitc polynomials. 

The Pq term will not make any contribution to first order 
because it is an odd function and thus, it gives a zero by 
(36) integrating over a Gaussian integral. On the other hand, 
the Pq term will make a nonzero contribution to first 
order. The contribution of the Pq term to first order is 
given by 



^-^OiGUP) 

and thus we get 



|L<^„|,.iL|^„> ,44) 



AE, 



(1) - ^^ ri) V/ dx e--^ (3 -67x2+ 7 V) 



O(GUP) 



2m 



(45) 

where 7 is equal to ^. 

By integrating, we get the shift of the energy to first order 
of perturbation as follows 



AE, 



(1) 



15 



fc2 2 2 

n uj ma 



or, equivalently, 



A^, 



(1) 



15, 



E, 



(0) 4 



(46) 



(47) 



We now compute the contribution of the p'^ term to sec- 
ond order of perturbation 



AM2) 



y- |<^/c \Vi I V«>P 



E^'' E. 



(0) 



(48) 



III. SIMPLE HARMONIC OSCILLATOR 

We now consider a particle of mass m. The Hamiltonian 
of the simple harmonic oscillator with the GUP-corrcctcd 
Hamiltonian assumes the form 



where 



m dx^ 



(49) 



In particular, we are interested in computing the shift in 
the ground state energy to second order 



„2 1 

2m 2 



m 



Po 



5q:" 



■Po • (40) 



^^ = E 



\<xbk\Vi\^o>\' 



k^n 



E 



(0) 



E. 



(0) 



(50) 



and for this reason we employ the following properties of 
the harmonic oscillator eigenfunctions 



< -ipjn \ X \ Ipn >-- 



and 







TO 7^ n ± 1 
^ , TO = n+l 
^ , TO=n- 1 



(51) 



measurement in the case of J/ip [15] is at the level of 
10~^. This sets the upper bound on ag to be 



ao < 10^ 



(58) 



It should be stressed that this bound is in fact consistent 
with that set by the electroweak scale. Therefore, it could 
signal a new and intermediate length scale between the 
electroweak and the Planck scale. 



< i>m \ X^ \ IpO > = ^ < 1pm \ X \ 'ipk >< Ipk \ X \ Ipl > 



k,l 



< %l)l \ X \ ll)Q > 



(52) 



which is nonvanishing for the {l,k,m) triplets: (1,0,1), 

(1,2,1), and (1,2,3). 

Thus, the ground state energy is shifted by 



(2) _ a^h'^ 



^K'' = 



771 •" 



^ |<^m|igT|V'0>P 



m^O 



E, 



(0) 



E^, 



(0) 



1/4 



Since the eigenfunction | -00 >= (■^) e "r^ , we 
have -^ \ ipo >= (37^2; - -f^x^) \ V'o >• By employing 
these into Eq.(50), we get 



AE, 



(2) 



m#0 



\<lPm I (3x-7x3) I V;o >|- 

p(0) p(0) 



.(54) 



Using Eqs.(51) and (52), the energy shift finally takes the 
form 



a4^) = -H,^to(4")^' 



(55) 



IV. THE LAMB SHIFT 

For the Hydrogen atom, V{f) = —k/r {k = e^/47reo = 
ahc, e — electronic charge). To first order, the perturbing 
Hamiltonian Hi, shifts the wavefunctions to [16] 



\lpnlm)l = \lpnlm) + ^ 



^n'Vm'\nlri 



{n' I' m' }^{nlm} 



En 



(0) 



i^l?^ 



\tpn'l'm') (59) 



where n, /, m have their usual significance, and 

^n'l'm'lnlm = \Wn'l'm'\tl\\lpnlm) ■ 

Using the expression p^ = 2m[Ho + k/r] [8], the perturb- 
ing Hamiltonian reads 



Hi = -(aVSm) 
So for GUP effect to a order, we have 



r fci 




\ k] 


Ho + - 




ffo + - 


i ^ \ 




L ^J 



(60) 



en'l'm'lnlm = (V'n'/'m' | ( ) PoPo I V'nim) • (61) 

It follows from the orthogonality of spherical harmonics 
that the above are nonvanishing if and only if /' = Z and 



or, equivalently. 



AE, 



(2) 



11 



E, 



(0) 



— -Oi'-mE'"' = - 



(0) 



11 



mhuj 



(56) 



It is noteworthy that there are some systems that can 
be represented by the Harmonic oscillator such as heavy 
meson systems like charmonium [14]. The charm mass is 
TO-c ~ 1.3 GeV/c^ and the binding energy lo of the system 
is roughly equal to the energy gap separating adjacent 
levels and is given hy huj Ri 0.3Gel^. The correction due 
to GUP can be calculated at the second order of a. Using 
Eqs.(47) and (56), we found the shift in energy is given 

by 



62001100 = 2iafi {■)p2ao\ 



Ho + - 



§-) I'Aioo/ 



(62) 



We utilize the following to calculate the shift in the 
energy: 

(i) the first term in the sum in Eq.(59) {n' = 2) domi- 
nates, since En = —Eq/h? { Eq ~ e'^ /Sireoao = k/2ao = 
13.6 eV , ao — ineoh'^/me'^ — 5.3 x 10^^^ metre , to = 
electron mass = 0.5 MeV/c^), 

(fi) ipnim{r) = Rnl{r)Yim{0,(j)), 

(iii) i?io = 2a(7^/^e-'~/°" , and 

i?20-(2ao)-'/M2-r/ao)e-'^/2-o, 

(iv)yoo(e,<^) = l/(V4^) . 



AE, 



(2) 



E, 



(0) 



^0 






2.7 X 10 



-''a?^ 



(57) 



Once again, assuming oq ^ 1 renders the correction too 
small to be measured. On the other hand, if such an 
assumption is not made, the current accuracy of precision 



Thus, we derive 



6200 1 100 



2ahk 

-« (^^200 I 

ao 



IV 



100/ 



. SV^ahk 
' 27al 



(63) 
(64) 



Therefore, the first order shift in the ground state 
wavefunction is given by (in the position representation) 



where k = ^j2mE/h? and fci = ^2m(Fo' - E)/K 
Considering solutions of the form '0<,> 



e™^, we get 



A^ioo(f) ^ ^^ok^1-<i)(^1 



100 V 

32V2ahk , ^^ 

= ^ „, o „ V'200(?') 



E, -E; 



■^200 (^) 



Slag Eo 

. 64V2ah 

81ao 



^^200 (7^) 



(65) 
(66) 



Next, we consider the Lamb shift for the n*'* level of the 
hydrogen atom [17] 



Ai?W-|4fln-)|^„/„(0)p . (67) 

Varying V'nim(O), the additional contribution due to GUP 
in proportion to its original value is given by 



AE 



n{GUP) _ „A|V'n;m(0)| 



Thus, for the ground state, we obtain 



(68) 






AE- 



(1) 



64/i ao 
SlaoMpic 

1.2 X 10^22 



ao 



(69) 



The above result may be interpreted in two ways. First, 
if one assumes ao ~ 1, then it predicts a nonzero, 
but virtually unmeasurable effect of GUP and thus of 
quantum gravity. On the other hand, if such an as- 
sumption is not made, the current accuracy of precision 
measurement of Lamb shift of about one part in 10^^ 
[8, 18], sets the following upper bound on ag: 



ao < 10 



10 



(70) 



It should be stressed that this bound is more stringent 
than the ones derived in previous examples [9], and is 
in fact consistent with that set by the electroweak scale. 
Therefore, it could signal a new and intermediate length 
scale between the electroweak and the Planck scale. 



V. POTENTIAL STEP 

Next, we study the one-dimensional potential step given 

by 



V'(x) = K 0(x) 



(71) 



where 9{x) is the usual step function. Assuming E < 
Vq, the Schrodinger equation to the left and right of the 
barrier are written, respectively, as 









2iah(fip< = 
2iah(f'ip> = 



(72) 
(73) 



m^ -f fc^ -f 2iahm^ == 
m^ - kl + 2iahm^ = 



(74) 
(75) 



with the following solution sets to leading order in a, 
each consisting of three values of m 



x<0 :m = {ik' ^ ^ik" , —-} (76) 

2ah 

x>0 : m == {fc'i, -fc'/, -^} (77) 



where 



k' = fc(l + kah), k" = k{l - kah) (78) 

k[ = ki{l-iahki), k'{ ^ ki{l + iahki) . (79) 

Therefore, the wavefunctions take the form 



Ae*'^'^ + Be-'''"'' + Ce^ , x<0 (80) 



ipy = De-''"'' +Ee^, 0<x 



(81) 



where we have omitted the left mover from i/'> . 

Now the boundary conditions at a; — consist of three 

equations (instead of the usual two) 

d"V<|o = rf>>|o, n = 0,l,2. (82) 

This leads to the following conditions: 

A + B + C = D + E (83) 

(84) 



i(k'A-k"B + —) = -k'{D + — 



k'^A + k"^B 



2ah' 
C _ 

{2ahf ^ {2ahf 



E 



2ah 
kfD . (85) 



Assuming C ^ E ^ ©(a^), we have the following solu- 
tions to leading order in a 



B _ ik' + k'{ 

A ~ ik" - k'{ ' 

D _ 2ik 

'a ^ ik" - k'{ ' 

E-C k'^ (ik" - k'l) + k"^ [ik' + k'{) + fcf {2ik) 



(86) 
(87) 



{2ah)^A 



ik" - k'{ 



It can be easily shown that the GUP-corrected time- 
dependent Schrodinger equation admits the following 
modified conserved current density, charge density and 
conservation law, respectively, [9] 



2rai \ dx dx 



ati^ /d2|^|2 dtpd^j* 

I 1 : ~ -^ 



m V dx^ 



dx dx 



I ;|2 dJ dp 

1^1 ' a^+^F-'^- 



(89) 
(90) 



The conserved current is given as 



m ^ 
2ah^k^ 



\B[ 



{\A? 



\B\ 



\Cf 



(91) 



m am 

The reflection and transmission coefficients are given by 



R 



B 

A 

ik' 



1 - 2ahk 
1 + 2ahk 



^k'{ 

(fc2. 



' 1 - 2ahk 
1 + 2ahk 



1 - 2ahk 



{kl + k'^Y{l-Aahk)l + 2ahk 



T 



\jj\2^-2k,x 



|^|2g-2feia; 



hk 



\A\^{l + 2ahk) 



0, 



R + T = 1. 



(92) 

(93) 

(94) 

(95) 



At this point we should note that GUP did not affect R 
and T up to 0{a). 



VI. POTENTIAL BARRIER 

In this section we apply the above formalism to an STM 
and show that in an optimistic scenario, the effect of the 
GUP-induced term may be measurable. In an STM, free 
electrons of energy E (close to the Fermi energy) from 
a metal tip at x = 0, tunnel quantum mechanically to a 
sample surface a small distance away at x = a. This gap 
(across which a bias voltage may be applied) is associated 
with a potential barrier of height Vq > E [19]. Thus 



V"{x) = V^' [e{x) - 9{x - a)] 



(96) 



where 9{x) is the usual step function. The wave functions 
for the three regions, namely, x < 0, < x < a, and 
X > a, are iI'1j4'2j and ips, respectively, and satisfy the 
GUP-corrected time-independent Schrodinger equation 

rf^V'1,3 + fc^V'1,3 + 2iaM^ipi^z = 
cPi'2 - kf'ip2 + 2iahcP'il)2 = 



where k = ^j2mE/h? and fci = yj2m{V^' - E)/h? . 
The solutions to the aforementioned equations to leading 
order in a are 



V'l 


= Ae^^'^ + Be^"^^"^ + pg*^/^"'! 


(97) 


V'2 


= Fe''!^ + Ge-'^"^ + Qe"/^"'^ 


(98) 


■fAs 


= Ce^^' ^ + i?e*^^^"'' 


(99) 



where k' — fc(l + ahk),k" = fc(l — ahk),k[ — fci(l — 
iahki),k'( = fci(l + iahki) and A,B,C,F,G,P,Q,R are 
constants of integration. In the above, we have omitted 
the left mover from ip^. Note the appearance of the new 
oscillatory terms with characteristic wavelengths ^ ah, 
due to the third order modification of the Schrodinger 
equation. The boundary conditions at x = 0,a are given 

by 



d"'ip2\x=a = 

If we assume that P 
following solutions 



d"^2U=o , n = 0,l,2 (100) 

rf"^3U=a , n = 0,l,2. (101) 

-^ Q r^ R ^ ^(a^), we get the 



C _ i{k'k'{ + k"k[ + k'k[ + k"k'{)e-'''' ''+''> 

A ^ f,{k'^+k'()a(j^, _^ jfc^)(fc// + ik'r^ _ (fc/ _ ik'l){k" - ik[y 

(102) 

(103) 
(104) 



(105) 



B 

A 

F 


kl + ik ik'a-k\a^ 
k'{ - ik" [ A 

[l+i^)e^k'a-k\a^ 


-1 


A 
G 


(l^ifc;)e*'='''+fe""§ 





A 



1 + 



From Eq.(89), it follows that the transmission coefficient 
of the STM, given by the ratio of the right moving cur- 
rents to the right and left of the barrier, namely, Jr and 
Jl, respectively, is to 0{a) 



T = — 
Jl 



r. 


2 


B 




- 2ahk 




A 




A 



(106) 



which gives using the solutions in Eqs.(102) and (103) 
the following final expression 



T = To [1 + 2ahk{l - T-^)] 



Tn = 



1&E{V^'-E) 



2kia 



K 



112 



(107) 
(108) 



where Tq is the standard STM transmission coefficient. 
The measured tunneling current is proportional to T 
(usually magnified by a factor Q), and using the following 
approximate (but realistic) values [19] 



■m^me= 0.5 MeV/c^ , E ; 



Vo" = 10 ey 



a= 10 



-10 



/o = 10-'' A , g = io^ 



we get 



/o To 



SI = gsi = 10 



-26 



A 



(109) 



where we have chosen oq = 1 and Tg = 10 ^, also a 
fairly typical value. Thus, for the GUP-induced excess 



current 51 to give the difference of the charge of just one 
electron, e ~ 10~^^ C, one would have to wait for a time 



^-Jl-'' 



(110) 



or, equivalently, about 4 months, which can perhaps be 
argued to be not that long. In fact, higher values of 
ao and a more accurate estimate will likely reduce this 
time, and conversely, current studies may already be able 
to put an upper bound on ao- 

What is perhaps more interesting is the following re- 
lation between the apparent barrier height <^a = Vq — E 
and the (logarithmic) rate of increase of current with the 
gap, which follows from Eq.(107) 



h 


din/ 


VSm 


da 



8m(kki) 



(111) 



Note the GUP-induced deviation from the usual linear 
^y^A vs \dhil/da\ curve. The exponential factor makes 
this particularly sensitive to changes in the tip-sample 
distance a, and hence amenable to observations. Any 
such observed deviation may signal the existence of GUP 
and, thus, in turn an underlying theory of quantum grav- 
ity. 



VII. CONCLUSIONS 

In this work we have investigated the consequences of 
quantum gravitational corrections to various quantum 
phenomena such as the Landau levels, simple harmonic 
oscillator, the Lamb shift, and the tunneling current in a 
scanning tunneling microscope and have found that the 
upper bounds on ao to be 10^"^, 10^^, and lO^'^ from the 
first three respectively. The first one gives a length scale 
bigger than electroweak length that is not right experi- 
mentally. It should be stressed that the last three bounds 
are more stringent than the ones derived in the previous 
study [9], and might be consistent with that set by the 
electroweak scale. Therefore, it could signal a new and in- 
termediate length scale between the electroweak and the 
Planck scale. On the other side, we have found that even 
if ao ~ 1, we still might measure quantum gravitational 
corrections in a scanning tunneling microscopic case as 
was shown in Eq. (110). This is in fact an improvement 
over the general conclusion of [9], where it was shown 
that quantum gravitational effects are virtually negligi- 
ble if the GUP parameter /3o ^ 1, and appears to be a 
new and interesting result. It would also be interesting 
to apply our formalism to other areas including cosmol- 
ogy, black hole physics and Hawking radiation, selection 
rules in quantum mechanics, statistical mechanical sys- 
tems etc. We hope to report on these in the future. 
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VIII. APPENDIX 

A. Proof for Eq. (1) 

Since black hole physics and string theory suggest a mod- 
ified Hcisenberg algebra (which is consistent with GUP) 
quadratic in the momenta (see e.g. Ref. [1] ) while DSR 
theories suggest one that is linear in the momenta (see 
e.g. Ref. [2] ), we try to incorporate both of the above, 
and start with the most general algebra with linear and 
quadratic terms 

[xi,pj] = ih{Sij + Sijaip + a2^-^ + jSiSi^p^ 

+ P2P^P,) ■ (112) 

Assuming that the coordinates commute among them- 
selves, as do the momenta, it follows from the Jacobi 
identity that 

- [[xi,Xj],pk] = [[xj,pk],Xi] + [[pk,Xi\,Xj] = . (113) 

Employing Eq.(112) and the commutator identities, and 
expanding the right hand side, we get (summation con- 
vention assumed) 

= [[xj,pk],Xi] + [[pk,Xi],Xj\ 

= ih{-ai5jk[xi,p] - a2[xi,pjPkP~^] - l3iSjk[xi,pipi] 
-/32[xi,PjPk]) - {i^ j) 

= ih[-ai5jk[xi,p] - (X2{[xi,pj]pkP^'^ +pj[xi,pk]p^^ 
+PjPk[xi,p^^]) - l3iSjk {[xi,pi]pi +pi[xi,pi]) 
~P2 {[xi,pj]pk +pj[x,,pk\)) - {i^ j) . (114) 

To simplify the right hand side of Eq.(114), we now eval- 
uate the following commutators 

(i) [xi,p] to 0{p) 



Note that 



,p2] = [xi,p-p] = [xi,p]p + p[xi,p] 



\x,,,p-\ = [Xi,p- p\ = [Xi,p\p + p[Xi,p\ (115) 

[Xi,PkPk] = [Xi,Pk]Pk +Pk[Xi,Pk] 

ih (Sik + aipSik + a2PiPkP^^) Pk + ifiPk {Sik 
+aip5ik + a2PiPkP^^) (to 0{p) using (112)) 
2ihp,[l + {ai + a2)p] . (116) 



Comparing (115) and (116), we get 

[xi,p] = ih {pip^^ + (ai + a2)pi) . 

(ii) [x„p-i] to 0{p) 



(117) 



Using 

= [xi,I] = [xi,p-p^^] = [xi,p]p^^ +p[xi,p^^] (118) 

it follows that 

[xi,p^^] = -Pi^[xi,p]p^'^ 

= -ih p^^ {piP^^ + (ai + a2)Pi) P^^ 

= -ih p,p-^ {1 + {ai + a2)p) . (119) 

Substituting (117) and (119) in (114) and simplifying, we 
get 

= [[xj,pk],Xi] + [[pk,Xi],Xj] 

= {ia^-a2)p-^ + ial+2|3^-(32))A,k^ (120) 

where Ajki — PiSjk — Pj5ik- Thus one must have ai = 
a2 = —a (with a > 0; The negative sign follows from 
Rcf. [3] of our paper), and /32 — 2/3i + a^. Since from 
dimensional grounds it follows that /3 ^ a^ , for simplicity, 
we assume /3i — a^. Hence /?2 == 3a^, and we get Eq.(l) 
of this paper, namely, 



[xi,pj] = ihi Sij-aipSij + -^-^ \ +a^{p^5ij + ipiPj) j .( 



B. Proof for Eq. (5) 



We would like to express the momentum pj in terms 
of the low energy momentum poj (such that [xi,poj] — 
ihSij). Since Eq.(121) is quadratic in pj, the latter can 
at most be a cubic function of the poi . We start with the 
most general form consistent with the index structure 



where a ^ a and b ^ a? . From Eq.(122) it follows that 

[xi,Pj] = [xi.PO] + apopoj + bplpQj] 

= zMy + a {[xi,p(i]poj + po[xi,po.j]) 

+ b{\xi,po]popoj+po[x^,po\poj+pl[x^,poj\) .(123) 



Next, we use the following four results to 0{a) and 
[xi,Poj\ — ifi in Eq.(123): 

(i) [xi,po] — "ifiPoiPo , which follows from Eq.(117) when 
ai — 0, as well from the corresponding Poisson bracket. 

(ii) Pj = Pqj{1 + apo) + C(a^) ^ Poj{^ + ap) [from 
Eq.(122)]. Therefore, poj ~ j|^ ~ (1 - ap)p, . 

(iii) po = {pojPoj)^ = ((1 ^ o^pfPjPj) = (1 - a.p)p . 

(iv) PoiPo^Poj = (1 - ap)Pj(l - ap)-^p-^{l - ap)pj = 
(1 - ap)piPjp-'^ . 



121) Thus, Eq. (113) yields 



[xt,pj] = ihSij +iah{p5ij +piPjP ^) 

+ ih{2b - a^)piPj + ih{b - a'^)p^5ij . (124) 
Comparing with Eq.(121), it follows that a — —a and 
b = IcP' . In other words 

Pi = Poj - apoPoj + 2a^PoPoj = P03 (l - apo + 2a^pl) (125) 



P'j = Poj + apoPoj + bplpoj 



(122) which is Eq.(5) in this paper. 
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